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Abstract 

The Wilsonian renormalization group (RG) method is applied to finite temperature 
systems for the study of non-perturbative methods in the field theory. We choose the 
O(N) linear sigma model as the first step. Under the local potential approximation, we 
solve the Wilsonian RG equation as a non-linear partial differential equation numerically. 
The evolution of the domain is taken into account using the naive cut and extrapolation 
procedure. Our procedure is shown to yield the correct solution obtained by the auxiliary 
field method in the large N limit. To introduce thermal effects, we consider two schemes. 
One in which the sum of the Matsubara frequencies are taken before the scale is introduced 
is found to give more physical results. We observe a second order phase transition in both 
the schemes. The critical exponents are calculated and are shown to agree with the results 
from lattice calculations. 

1 Introduction 

Non-perturbative aspects of the quantum field theory are very interesting and important. Study 
of low-energy hadron properties in QCD is an example that non-perturbative approaches are 
indispensable. We expect a phase transition at finite temperature which can be understood 
only non-perturbatively. Various non-perturbative approaches have been developed, such as 
1 E-mail address: umekawa@th.phys.titech.ac.jp 



1 



the super-daisy summation of diagrams, the instanton approach, 1/N expansion, e expan- 
sion and so on. The program, however, has not been accomplished yet. Under this circum- 
stance, the method of the Wilsonian renormalization group (RG) equation]!], ^ |3j receives 
much attention[|], |[ ||, [7|, ||, || . The idea is to include the quantum effects gradually from high 
energy in order to obtain a low-energy (Wilsonian) effective action. This RG equation is a non- 
linear evolution equation with respect to the scale A in the functional theory space. While A is 
lowered, various irrelevant operators are reproduced and finally we obtain the effective action 
in which all quantum effects are included. This process can be achieved non-perturbatively. 
The concept of the Wilsonian RG equation is common to various quantum field theories in the 
sense that it is model independent in contrast to the instanton or 1/N expansion approaches. 
A systematic approximation is also possible in the form of enlargement of the functional theory 
space step by step. 

In spite of these advantages, there are still technical difficulties in application. It is not 
easy to obtain the physical quantities while it is very useful to see the phase structure and 
the critical phenomena. This difficulty is especially serious in the symmetry broken phase. 
In order to treat the bottom flatened effective potential, which is naturally required in the 
quantum system, one must prepare a suitable functional theory space. There is a dilemma that 
the expected potential is non-analytic while the RG equation is a partial differential equation 
which assumes the smoothness of the solution. 

The aim of this paper is to challenge this dilemma. We here propose to solve the RG 
equation as a difference equation in the explicit scheme. In the symmetry broken phase, we set 
the domain in which the solutions is continuous, and the domain size is modified at each step 
of integration. This approach may sound too naive but we demonstrate the availability of such 
numerical treatment of the Wilsonian RG in physical applications. 

For this purpose we concentrate on the O(N) linear sigma model. This model is often 
studied in other non-perturbative methods and its physical application is also very interesting. 
It can describe not only the spontaneous symmetry breaking but also the symmetry restoration 
above a critical temperature. Especially, the 0(N = 4) symmetric model is known to be 
equivalent to the chiral SUl{%) X SUr{2) effective theory of QCD, and is widely expected to 
describe the chiral phase transition in QCD. 



2 



In Sect. 2, we present our formulation of the Wilsonian RG equation and the numerical 
method for solving the equation in the local potential approximation. We demonstrate that 
our method works very well for the large N limit, for which the auxiliary field method yields 
the correct answer. We apply our method for N = 4 and show that the symmetry breaking is 
well represented in our method. 

In Sect. 3, we apply our method to finite temperature O(N) model. We employ the imaginary 
time formalism. When the sum of the Matsubara frequencies is considered, we realize two 
possible schemes of treating the scale dependence. One is to cut-off the sum according to the 
scale change and the other is to take the sum before the scale is introduced. These two schemes 
should give the same results provided that the initial scale A is infinite. However, in practice 
we choose a finite Ao and the two schemes give different results. We examine both of them and 
in the end show that the second scheme is more appropriate. 

The results for the O(N) sigma model are presented. We obtain the second order phase 
transition to the chiral restored phase in the above two schemes. It is shown that the four 
dimensional system near the critical temperature in the first scheme reduces to a three dimen- 
sional system at zero temperature. This indicates the first scheme corresponds to taking the 
high temperature limit even near the critical temperature. On the other hand, in the second 
scheme we obtain more physical results. The critical exponents are computed and are shown 
to agree very well to the results in the lattice calculation. 

The conclusion is given in Sect. 4. 



2 Wilsonian RG Equation 

The lagrangian of the O(N) linear sigma model we consider here is given by 

in the Euclidean notation /i = 1,2,3 and 4. <fi denotes the A" component column vector 

<P ■= (01,02, ■ • ■ AnY- 

The Wilsonian RG equation with a sharp cut-off scheme is called Wegner-Houghton (WH) 
equation. The WH equation in the O(N) linear sigma model is derived in Ref. JIJ and is given 
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by 



dr\ = Af J / 1 ^r A \ ar A / £ 2 r A \ 1 <yr A 

dA 2 4i=a P l 1 n \A 2 6(/> p 6<t>- p ) % \5(j> p 5(j>. p ) 50 V l_) 



in .D dimensional space-time. r A [0] denotes a scale dependent effective action at the scale A 
above which the quantum fluctuation is included. Thus at A = the scale dependent effective 
action agrees with the conventional effective action r [</>]. The integral of r.h.s. is the shell 
(\p\ = A) integral 

The first and the second terms in r.h.s. of Eq.@ give ring and dumbbell diagrams respectively. 
For a given initial condition r Ao [0] at the scale A , the effective action T A [0] at any scale A is 
obtained by solving the WH equation. As the initial condition for Eq.(0) we assume that the 
effective action T A [0] reduces at high A to the classical action S [</>]. 

In practice, to solve the equation (§), suitable approximation is necessary. We employ the 
local potential approximation (LPA) which assume the functional space as 

r A [0]= J ^x^-id^f + v^y (4) 

Under this approximation, the dumbbell diagram does not contribute and the WH equation 
reads 

where A := A exp(— t), x := \/<P, 

n -D/22l-D 

Ad ■= w (6) 

and V t (x) is a scale dependent effective potential. 

We solve this equation numerically as a difference equation in the explicit scheme 

^ = {(TV - 1) in f 1 + 

At 2 \ v J \ tijxi 2Ax J 

+ 1 H 1 + A| J I (7) 

where tj := jAt, Aj := A exp(— tj), Xj := iAx and V? := V tj (xi) with j = 0, •••,j max and 
i — 0, • - • , i max . This explicit scheme suggests us that a large j max is necessary for a large z max . 



Furthermore there is a subtle problem in solving Eq.(^), that is , we have to specify the domain 
D of x. In the broken phase, the field variable can not take all the values because the field 
variable is not an analytic function of the source J around the origin = 0. Numerically the 
arguments of logarithm become zero or negative for x with \x\ < 3 a at some tj. In such cases, 
we obtain the effective potential at only by the information of the region x with \x\ > a. 
Thus, we consider the domain T> = {x | \x\ > a}. On the other hand, when the arguments of 
the logarithm grows up from zero at the end point x = a of the domain D at tj, we extrapolate 
the arguments as a smooth function of x using vicinal four points by 

f{Xm-l) ~ f{%m) - f{x m )Ax 

~ f{xm) ~ {f'{x m+1 ) - f"(x m+1 )Ax) Ax 

f( v ( f{Xm+2) ~ f{x m ) f(x m+2 ) ~ 2f(x m+ l) + f(x m ) \ ^ , g x 

J[Xm> \ 2 Ax ~ Ax ) X U 

where x m is the nearest neighbor to a and we enlarge the domain T> at t J+1 . The numerical 
result seems to be almost independent of the extrapolation procedure if we take sufficiently 
large i max (and j ma x) • 

We choose the mass unit as A so that A = 1, concentrate on D = 4, and set \i 2 = /i 2 /Ag = 
—0.6 and X/8N = 4 so as to bring out the spontaneous symmetry breaking. We employ the 
numerical parameters z max = 400 and j max = 10 5 . 

First we consider the large N limit, where the LPA WH equation (|^) reduces to 

dt 2 \ A 2 x dx J K ' 

Eq.(|9D is solved for the initial condition 

V t=0 (x) = ^X 2 + A x 4_ (1Q) 

The result is shown in Fig.|I]. For t = 0.0 ~ 0.5, the effective potential becomes gradually flat 
as the quantum correction is included. This resembles the lattice result in Ref . |T0| . For t > 0.5, 



however, the effective potential becomes complex from the inside and the Domain develops. As 
stated above, we do not use the values of such complex region in computing the next step. It is 
plausible that this region should be perfectly flattened as is indicated by the study at t — > oo in 
Ref. pi- It should be noted that the point of the potential minimum, which gives the vacuum 



5 




Figure 1: The solutions of the Wilsonian RG equation for the large N limit. The curves 
correspond to t — 0.0, 0.2, 0.5, 1.0, 9.9 from bottom to top respectively. 

expectation value of 0, and the end point of the domain become close to each other. Finally, at 
t — ► oo, these two points coincide and the effective potential defined conventionally is obtained. 
It is well-known that the exact effective potential can be obtained using the auxiliary field 



method in the large N limit [pi] , [12|j . As a result, the scale dependent form of the effective 
potential is given by 

V t (x) = V t (x,s) = -x 2 s-—s 2 + ^s + — ——\n(k 2 + s) (11) 



2 2A A 2 J a (2tt) 4 
where the auxiliary field s is a function of x and is given by a solution of the saddle point 

condition 

= 0. (12) 

OS 

It is easy to see that Eq.flTTD is the solution of Eq.([J) with the initial condition (|1C|) |T3| . We find 
that our numerical solutions in Fig.0 agree almost completely to the analytic solution Eq.(ll 



In plotting Fig.0, we only take the region of <ft where V t (x) is real. In fact, for each t there 
is a critical value of x c (t) = \4> c (t)\ below which the solution of Eq. fljjD becomes complex. In 
the present case, 

Xc (t = oo) = 2N nun | S -^- - ( A 2 - s In ^^-) \ (13) 
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Figure 2: Solution of the Wilsonian RG equation for N = 4. The curves correspond to 
t = 0.0, 0.5, 10 from bottom to top respectively. 

This phenomenon occurs where the O(N) symmetry is spontaneously broken and the system 
generates a non-zero vacuum expectation value of 0. The magnitude of the vacuum condensate 
of 4> is given by the critical value at t — ► oo or A — > 0. 

Next we consider the finite N case. The 0(4) result is shown in Fig.^. Again the complex 
region appears. One sees that the qualitative behavior of the effective potential is quite similar 
to that for the large N limit and that the calculation supports the spontaneous symmetry 
breaking. The evolution of the domain by the scale is shown in Fig.[| We see that the curves 
are monotonically increasing and saturate at about t ~ 5 (A ~ 0.007). 




3 Finite Temperature 

The introduction of the thermal effects is discussed in many text books |14]]. We use the 
imaginary time formalism in which the loop momentum integration is replaced by the sum of 
the Matsubara frequencies as 

j^o^T^j^, Po = 2n«T (14) 
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Figure 3: The end points a of the domain v.s. the scale variable t. The curves correspond to 
the temperature T = 0.0, 0.1, 0.15 from top to bottom respectively. 

where T denotes the temperature. The system is now controlled by the new dimensional 
parameter T. There are some candidates how to extend the evolution equation in order to 
include the effect of T. We consider the following two schemes within the local potential 
approximation. 

Scheme I 

Firstly we consider the 4-dimensional spherical cut-off scheme. For T = 0, the formula 

A 1 f A 1 f A J D-l f dQ P 



5AJa-8A< p <a(2tv) d SAJasa J (2n) D 

J (2ti) d J\ p \=a y K J 

is used in the derivation of the WH equation (Q). If a function f(p 2 ) depends only on p 2 , then 
Eq.flT5"D reduces to 

/ . . i^f(P 2 ) A D ^d/(A 2 ) (16) 



8 A Ja-sa<p<a (2tt) d ' 

which is used as the prefactor in r.h.s. of Eq.([|). Even at finite temperature, the integrand is 

still f(p 2 ) within LPA and therefore the shell integration can be performed as 

11 r dP-^y 2 

8kp^h-8K<\p\<A {2ti) d - iJ{P ' 



OA p JA-.5A<- v /(27rn//3) 2 +pp<A 



A 1 _ rv /A2 -( 2 ™// 3 ) 2 

|„|<aV/(2 7 t) V( a -« 5A ) 2 -( 2 ™/^) 





^° E l-(^) ) /(A 2 ) (17) 

^ |n|<A/3/(27r) \ V ^ 

where /3 := 1/T. Comparing Eq.(|16|) and Eq. fllTD , we conclude that the thermal effect is taken 
into account by replacing the prefactor of r.h.s. of Eq.(|^), 

A D K°^A D ^— £ 1- Kjl) • (18) 

Finally, we obtain the WH equation in the scheme I, 
dV t (x) _ A D ^k D - 1 _ / /27rn\\^ 

Here it should be noted that the thermal loop contribution is summed up only for \n\ < KqI3/2tt 
in the scheme I. For Ao — > oo, this scheme provides the full thermal loop correction, while for 
relatively small A , the summation is limited to low \n\ especially at high temperature. Indeed, 
we will see later that the limited summation gives a fluctuation in the effective potential. In 
order to avoid this difficulty, we next consider the second scheme. 

Scheme II 

We consider the scheme in which the cut-off is introduced after the separation of the quantum 
and thermal effects. In this scheme, therefore, we perform the summation of all Matsubara 
frequencies. 

For simplicity, we focus only on LPA. Suppose that the effective potential Vt(4>) is given at 
the scale t. Then the effective potential K+At(0) is given by 



= W) + ^E/^r^ 



2/5 (2tt) 



, 2 + fW) 



+ (2-loop and higher loop contributions) (20) 
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where the integral is performed between t and t + At. To specify the cut-off scheme we use the 
formula 

+ |/ ^p^r ln(l - exp(-3VF+IF)) 

= / dA\— ln(A 2 + M 2 ) + -5zL ln(l-exp(-/5v / A 2 + M 2 )) (21) 

ja=o 12 p J 

where (3 = 1/T and M 2 is an arbitrary constant. In r.h.s. the quantum and thermal effects 
are separated. Thus we introduce a sharp cut-off of A in Eq.(pl|). Similar procedure can be 
performed for the 2-loop and higher loop contributions. The cut-off is introduced so that in the 
quantum correction part the radial mode in the 4-dimensional integration is cut-off, while in 
the thermal correction part radial mode in the 3-dimensional integration which is left after the 
summation or the integration with respect to ko is cut-off. Taking the limit At — > in Eq . (|20|) , 
the 2-loop and higher loop contributions vanish. Therefore we conclude that the thermal effect 
is included by adding 




+ In [ 1 - exp A 2 + ^ | | \ (22) 

to r.h.s. of Eq.(||). 

Finally we obtain the WH equation in the scheme II. 



+ ln(l-exp(-^A 2 + ^] ] }. (23) 

It is noted that neglecting the quantum effects (the first term in r.h.s.), this evolution equation 
coincides with the thermal RG equation in Ref . |15j , which is derived in the real time formalism 



T6[], up to an additional step function. 



10 



In these two schemes, both Eqs.(|T^) and (p3| ) reduce to the zero temperature formula Eq.(|5|) 
in the limit T — > and they reduce in the limit T — > oo to the effective D — 1 dimensional form 

The difference between the two schemes come from the cut-off procedure. Taking A — > oo, 
their effective potentials at A = coincide. In practice, however, we neglect the thermal effect 
above Ao in the first scheme and assume that Ta equals to the zero temperature action. Thus 
the two schemes give different results for finite A . 

In the scheme I, the numerical results are as follows. Fig.f| shows the vacuum expectation 
value of as a function of the temperature T in the large N limit of Eq . (|T9|) . Fig.|5| shows the 
same quantity calculated by using the auxiliary field method. As these two curves agree well, 
we conclude that our numerical approach is appropriate. The result for A^ = 4 is given in Fig.^|. 
In these graphs, one sees unnatural steps. They come from the summation in Eq.(|l^). The 
summation is limited to a few terms when T becomes comparable to A /2it 150 [MeV] for 
A = 1 [GeV]). As the vacuum expectation value is sensitive to the number of the summation 
in the initial condition, there appear unnatural steps seen in the results. Since the critical 
temperature T c = 0.185 shown in Fig|4| is larger than Ao/27r, the phase transition in this 
scheme I is controlled by the high temperature formula Eq . (|24]) . At low temperature, we do 
not see such steps and the results agree very well with those obtained in our second scheme. 

In Fig.[7], we plot log(0) as a function of log(T c — T), where T c denotes the critical temper- 
ature. After fitting the curve to the linear function using the least square method, we obtain 
the critical exponent (3 = 0.448, defined by 

<f)(x\T c -Tf. (25) 

We also obtain the critical exponent 5 = 3.80 from Fig||, where 6 is defined by 

V oc \<f)\ 5+1 (26) 

at T = T c . In the same way we obtain the critical exponent 7 = 1.43 from Fig|| where 7 is 
defined by 
d 2 V 



S2 



oc |T C - T| 7 . (27) 
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Figure 4: Temperature dependence of the condensation as a function of temperature T in 
the scheme I for the large N limit. 

These results are summarized in Table [I]. From Eqs.(|25|)-(|27D, we obtain a scaling relation 
among the exponents, 



7 = 0(6-1). 



(28) 



Our numerical results yield (3(8 — l)/j = 0.877, which is not far from 1. 

Numerical results of the scheme II are shown in Figs.|lO|-|T2| and Table [l]. In this scheme 
the summation of the infinite number of the Matsubara frequencies is performed. Therefore 
in the scheme II the unnatural steps seen in the scheme I do not appear. These curves again 
show the second order chiral phase transition. We estimate the critical exponent (3 = 0.384 



from Fig.|13|. This value agrees with the value 0.384 obtained by the Monte-Cairo calculation 
in Ref. [TJ§, where the O(N) Heisenberg model in three dimension is used. We also obtain the 



critical exponent 8 = 4.70 from Fig. 14. The value given in Ref. [17] is 8 = 4.85. In the same 



way, we obtain the critical exponent 7 = 1.38 from Fig.|15|. The value in Ref.[I7| is 7 = 1.48 
which is also close to ours. These critical exponents satisfy the scaling relation very well, giving 
/3(<J-l)/ 7 = 1.03. 
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Figure 7: The critical exponent (3 for the result of the Wilsonian RG method for N = 4 in the 
scheme I. 
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Table 1: Critical Exponents 
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Figure 10: The temperature dependence of the condensation in the scheme II for the large 
N limit. 
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Figure 11: The same as Fig. 10 in the auxiliary field method. 
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Figure 12: The temperature dependence of the condensation in the scheme II for N = A. 
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Figure 14: The critical exponent 8 for N = 4 in the scheme II. 




4 Conclusion 



We have applied the Wilsonian RG method as the non-perturbative calculation method in the 
field theory to the O(N) linear sigma model and calculated the low energy effective potential 
in the local potential approximation. We have introduced a numerical recipe to define the 
domain for the effective potential and have tested the method in the large N limit. The results 
coincide with the solutions obtained by the auxiliary field method very well. This means that 
our numerical treatment of the Wilsonian RG equation is sufficiently good in this situation. 
We have also shown that the N = 4 effective potential behaves similarly as those in the large 
N limit. 

We have applied our approach to the finite temperature system to see the chiral restoration 
phase transition. We have employed the imaginary time formulation using the sum of the 
Matsubara frequencies. When we take the sum of the Matsubara frequencies, we consider 
two schemes. In the scheme I, we use the 4-dimensional spherical cut-off. Then the sum 
of the Matsubara frequencies is also limited accordingly. In the scheme II, we separate the 
quantum and thermal effects, before the momentum integral is cut-off. These two schemes 
should agree with each other if the initial Ao is infinite. In the scheme I, however, we assume 
that the thermal effect above Ao is negligible and T/y equals to the zero temperature action. 
Thus the two schemes give different results for finite A . Both of them show the second order 
phase transition. But the results of the scheme I have unnatural steps. These steps arise in 
the course of transfer from the four dimensional system to the three dimensional one at high 
temperature. On the other hand we have found that the scheme II gives more natural and 
physical results. Especially, the critical exponents in the scheme II are very close to that of 
Monte-Cairo calculation. 

Our aim has been to check whether our numerical method can be reliably applied to the 
Wilsonian renormalization group equation. We have demonstrated that it works very well in 
the O(N) linear sigma model. It has been shown also that the method can be applied to 
finite temperature system. This is the first step of our approach to more complicated systems 
including fermions and gauge fields. The Wilsonian RG method provides us with a powerful 
tool to study non-perturbative features of low energy effective theories. As a next step, for 
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instance, this approach may be able to judge whether certain effective model is consistent with 
QCD by comparing the solution of the Wilsonian RG equation for the effective model and the 
result of the lattice QCD. 
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